The interplay between double exchange, super- exchange, 
and Lifshitz localization in doped manganites 
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Considering the disorder caused in manganites by the substitution Mn^ Fe or Ga, we accomplish 
a systematic study of doped manganites begun in previous papers. To this end, a disordered model 
is formulated and solved using the Variational Mean Field technique. The subtle interplay between 
double exchange, super-exchange, and disorder causes similar effects on the dependence of Tc on 
the percentage of Mn substitution in the cases considered. Yet, in La2/3Cai/3Mni_yGay03 our 
results suggest a quantum critical point (QCP) for y « 0.1 — 0.2, associated to the localization of 
the electronic states of the conduction band. In the case of LaxCaajMni-^FeyOa (with x = 1/3, 3/8) 
no such QCP is expected. 

PACS numbers: 75.30.Vn, 71.30.+h, 75.50.Lk, 75.10.-b. 



I. INTRODUCTION 

A continuous phase transition at zero temperature de- 
fines a Quantum Critical Point (QCP), with a diverg- 
ing length-scale for quant um-pfkiet nations and a diverg- 
ing time scale for dynamics JJlo'cl The Renormalization 
GroupQ can be extended to the study of QCP,H which 
has become a popular concept among theoreticians, given 
the enhanced predictive power of the calculations. For 
instance, the experimental finding of a QCP has been 
postulated to uaderlie the physics-,p|f high-temperature 
superconductorsO and spin-ladders,EJIa it has been used to 
predict the scaling behavior for meta-magnetic quantum 
criticality in metals,LJ orJ^o study the effects of disorder 
in quantum spin-chains. □ The localization transition, of 
non interacting electrons in the presence of disorders is 
another type of continuous phase transition at zero tem- 
perature. The experimental finding of a QCP compatible 
with the theoreticabsxpectations for the Anderson metal- 
insulator transitionO would be very interesting. Note that 
the Anderson transition is not the cause of the metal- 
insulator transition in Lai_^CaxMn03jl3 a well known 
colossal magnetoresitance (CMR) manganite.til 

In this paper, we want to address the behavior of typ- 
ical CMR manganites when doped at the Mn sites 



Lai_:rAE:rMni_;^TRj^03 , 



(1) 



where AE = Ca, Sr and TR = Fe, Al, Ga, in the range 
0.3 < X < 0.5. Some expe rimc j iital . w ork has already 
been done for this materialsOEjEdiijIiB but, from the 
present analysis we believe that several interesting fea- 
tures have not still been found experimentally. Specifi- 
cally, we will show that, from the knowledge gained in the 
study of the phase diagrams of (Lai_^RE^)i_2;Sr:rMn03 
and (Lai_2RE2)i_a;Ca2;Mn03, where RE is a trivalent 



rare earth, and from our model calculations, one can con- 
clude that: 

i) When the y = Q material shows a first order tran- 
sition from the paramagnetic state (PM) tO|-the ferro- 
magnetic (FM) one, hke in La2/3Cai/3Mn03liZl, a QCP 
should appear close to 15% Ga substitution {x — 1/3 and. 
y = 0.15). Unfortunately previous experimental workO 
did not go bevond 10% of Ga, but more experiments 
will be done.Ea In this case, it is interesting that the lo- 
calization transition determines the magnetic behavior 
and then, through the corresponding magnetic critical 
exponents, it may be possible to investigate the physical 
properties of such QCP experimentally. In this context, 
note that disorder added to a system with a low T first 
order transition separating ordered phases (see Sec_|l^ 
and Fig. ^ produces quantum-critical like featuresO 

ii) When the y = {) material shows a second-order 
PM — FM phase transition, as in La2/3Sri/3Mn03, the 
low temperature phase will be ferromagnetic until values 
of y as large as w 0.4, and (if Mn is substituted with Ga) a 
FMmetal — FMinsulator transition at finite temperature 
should take place (see also Ref. 20). 

The layout of the rest of this paper is as follows. In 
the next section, we describe the experimentally known 
effects of Mn substitution and we will discuss the sim- 
plest model that can account for it. We will first consider 
the relatively simpler case of Fe Mn substitution, and 
we will then address the case of Gallium. From now 
on, we will not be concerned with the Al case, because 
this substitution causes the lattice to loose oxygen for 
y K, 0.1, and, a significant amount of oxygen vacancies 
are present .E3 This phenomenon can easily be taken into 
account with a proper modification of the present model, 
but this is left for future work. In section III we present 



our results, while the conclusions, comments, and predic- 
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tions are left for the last section. 



II. MODEL 



Fe substitution 



In previous works jBlnilsi'H we have considered the 
simplest model that describes the magnetic properties 
of (Lai_2RE^)i_2,AE^Mn03 and we concluded that, for 
0.3 < a; < 0.5 such model is the single orbiijal dou- 
ble exchange model (DEM) on a cubic latticec3 with 
isotropic hopping and super-exchange antiferromagnetic 
first-neighbors interaction: 



7^ = ^^(s.,5,K^ 



(2) 



where Ci corresponds to a Mn eg orbital while S stands 
for a unit vector oriented parallel to the Mn'it. {S — 3/2) 
core spins, which we assume to be classical.cJ The func- 



tion t{Si,Sj) — ^[cos^cos-^ -I- sin Y sin ^e 



stands for the overlap of two spin 1/2 spinors oriented 
along the directions defined by Si and Sj, whose po- 
lar and azimuthal angles are denoted by 9 and ip, re- 
spectively. The hopping integral, is approximately 
0.16 eVJli As far as the magnetic interactions are con- 
cerned, one can trade the effects cij3honons,LZI by tuning 
the superexchange constant Jaf-" The phase diagram 
of the model (Q) has been extensively studied by., means 
of the variational mean-field (VMF) teckaiqueEIIc3 and 
by the Hybrid Monte Carlo method.tJO These stud- 
ies have shown that the parameter Jaf is rather con- 
strained. On one hand, if one wants to reproduce the, 
first-order character of the PM-FM phase transition,EZl 
one should have Jaf larger than 0.06t.c2l But, on the 
other hand, the system is FM at low temperatures,only 
up to Jaf ~ O.OSi, for a hole concentration of 1/3.0 The 
conclusion is twofold: i) Jaf is not a very tunable param- 
eter and, ii) La2/3Cai/3Mn03 is a rather critical system 
where small perturbations can produce large changes in 
physical properties. 

Therefore, it is a challenge to this model to give 
account of the experimental situation of |-tl|ie,pjd(mblj3 
series (0) in the range 0.3 < a; < 0.5ll3Oll3By, 
without stretrhjug out the values of Jaf already 
determined H'EBo 

Through the whole series (|l|), only the Mn eg (up) 
band is electronically active, with the electron hopping 
between Mn^+ and Mn'^^. In the Fe case, the Fe eg 
(up) band is completely filled and electron hopping from 
Mn^+ to Fe'^"'" is forbidden. The situation is similar when 
TR = AP+ or Ga'^^, in which cases the d-band is full. 
This means that the fraction of Mn''^+ with respect to the 
amount of Mn^+-|-Mn*+ is incresed by a factor (1 — y)^^ 
respect to the case ?; = 0. Figure 2 of Ref. pfl shows 
that, as far as the the following argument is concerned, 
the Curie temperature, Tq does not depend much on the 



hole density around 40% of carriers. This will allow us 
to extrapolate the results for TR = Fe from our previ- 
ous computations. In fact, the similar ionic radii of Fe^"*" 
and Mn'^^ means that lattice distortion effects may be 
ignored. Thus, let us explore the hypothesis that the 
substitution of Mn for Fe affects only the antiferromag- 
netic interaction between the localized spins, besides the 
irrelevant in this case (1 — y)^^ factor in the effective 
number of holes. This change is due to the arising of 
Mn-Mn, Mn-Fe, and Fe-Fe couplings, with probability 
(1 — y)^, 2y{l — y) and y^ respectively. Given that the 
core spins are 3/2 for Mn and 5/2 for Fe, the effective 
value of the super-exchange constant is 



Jeff — Jaf 



(l-y)2 + ^2y(l-y) + |y2 



(3) 



Thus, we can easily extrapolate the results of Fig. 3 of 
Ref. |l| for a; = 3 /jS, to the case of Fe substitution (we 
took Jaf/^ ~ 0.07EJ), where we get Jeff — 1.17Jaf for 
y = 0.12, what implies that Tc(?/ = 0.12)/rc(2; = 0) = 
0.4 (see Fig. 3 of Ref. ^l|), in complete agreement both 
with the experimental data of Ref. |l^ as well as 'pith the 
phenomenological analysis of A. Tzavellas et aZ..E3 

On the other hand, when TR in ([^) is Al or Ga, lattice 
distortion effects enhance the relevance of the fact that 
on the Mn sites the charge is a mixture of +3 and +4. In 
fact, as we will see, electrostatic effects yielding diagonal 
disorder in our model turn out to be cru cial not only to 
understand the known experimental facts,t3t3ll3 but also 
to produce a new and very interesting physical situation 
for values of y not studied in Ref. O. 



B. Ga substitution 



The effects of Ga doping in the colossal magnetore- 
sistance materials La2/3Ca]^/3MnO, 



2/3'-^'l'l/3^ 

studied by preparing tha^series La2/3Ca2/3Mni_yVjaj,vj3 



has been recently 
^Gaj,03 

{y — 0, 0.02, 0.05, 0.1),I13 and a strong decrease of Tc 
with growing y was reported. However, the doping frac- 
tion never went beyond 10%, while we believe that a 
QCP will appear for y « O.L^O.2. No structural change 
upon doping was observedJ13 nor it is to be expected 
for larger Ga fraction, because it does not occur in the 
LaGa^Mni-yOs series. The model used here does not 
take into account the orbital degeneracy in the eg band. 
Thus, strictly speaking it ceases to be valid when the 
ratio Mn^+/Mn'^+ approaches unity, which, in this case, 
corresponds to y w 1/3. Nevertheless we use the results 
for y up to 0.4 in order to describe better the change in 
physical properties due to the transition. 

The Ga substitution disorders the system (details are 
given below). In order to describe specifically a disorder 
effect, we will need a disordered modelBB The standard 
way of constructing a disordered model is to choose first 
a simple model able to describe the ordered situation (i.e. 
the system with no Gallium), and then to mimic the in- 
clusion of impurities in the real system by the random 
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modification of some terms in the Hamiltonian that re- 
flect the microscopic effects of the impurity. As said in 
the previous subsection our reference ordered model, is 
the single orbital Double-Exchange Model (DEM) on the 
cubic lattice of Eq. (H) . 

With an ordered model in our hands, we can discuss 
the effects of Gallium doping. Given that the sizes of Mn 
and Ga are not very different, one expects small changes 
in the Mn-O-Mn angles, and therefore, on the values 
of the couplings t and Jaf- To keep the model sim- 
ple we shall assume that the parameter Jaf and t are 
not modified at all by Ga substitution. Since the series 
LaGaj^Mni„j^03 exist for all y, we expect Ga to be in a 
Ga^+ oxidation state, which is a fiUed-shell configuration. 
This means that the corresponding eg orbital will not be 
available for a hole in a neighboring Mn to hop, which 
reduces the Double-Exchange (DE) ferromagnetic inter- 
action. We thus encounter a quantum-percolation prob- 
lem (see e.g. Ref. ^0|). However, quantum-percolation 
studies have shown that the quantum threshold is not 
far from the classical one (that would be reached at 69% 
Ga fraction!), and this effect is not expected to be crucial. 
Moreover, given the fiUed-shell configuration of Ga, the 
Mn'^+ core spin 5'; will disappear from Ga-sites, which 
reduces also the AFM interaction that competes with the 
DE ferromagnetic one. In the case of Ga, the effective 
AF couphng of Eq. (||) is Jcs — — vYJaf which, for 
instance, for y = 0.12 turns out to be 0.77 Jaf- This is 
really a huge reduction (see Fig. 3 of Ref. 
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and Fig. 3 

of Ref. 24), and it is not obvious what the effect of Ga 
substitution on the Curie temperature would be. Fortu- 
nately, the effects of Ga substitution described so far cm, 
be straightforwardly studied with the VMF method,E3 
and it turns out that the reduction on the antiferromag- 
netic interaction has an stronger effect than the reduction 
of the DE mechanism, and Tq rises in the model when 
the Ga fraction is increased—Since this is in plain con- 
tradiction with experimentsjlj it is clear that our micro- 
scopic description of the Ga substitution is still incom- 
plete. This suggests that an on-site electrostatic pertur- 
bation is induced by Ga ions in Eq. (||) . This electrostatic 
term can be justified by the fact that the average charge 
on Ga sites is always -1-3 while on Mn sites it is a mixture 
of -f3 and -1-4. Therefore holes cannot hop onto Ga sites, 
but they are attracted to them. 

We expect this phenomenon to be more relevant in the 
case of Ga than in the case of Fe. In fact, if in a cubic 
lattice of Mn ions we substitute a fraction y of Mn by 
TR randomly, the electrostatic potential that appears in 
the Mn sites is proportional to the number of TR^+ that 
are neighbors of the considered Mn ion. Let us call e the 
electrostatic potential felt by a hole at a Mn site with 
only a neighboring TR. In the Fe case, e is the same on 
each Mn site because lattice distortion effects can be ig- 
nored, ought to the similar ionic radii of Fe'^"'" and Mn^+. 
However, in the Ga case (and even more in the Al case), 
the ionic radii difference implies that to the randomness 
in the distribution is added the randomness originated 



by the lattice distortion which, among other phenomena, 
makes the potential randomly dependent on the Mn po- 
sition. 

To avoid an overly complicated model, we incorporate 
the second disorder to the first one (of Lifshitz type) and 
we work with an effective e, the same for all Mn sites. 
Gonsequently, the e is larger for Ga and Al than feu: 
Fe; this is precisely what the experimental data requireE3 
since, as we have said, without a localization mechanism 
that destroys the double exchange, Tc would rise with y. 

As we will see, a moderate disorder (from our analysis, 
e < 2t) modifies only slightly the decrease of Tc with y. 
For this reason we have neglected it in the case of Fe. 

In summary, our model for La2/3Gai/3Mni_j,Gay03 is 
as follows: every Ga atom modifies the Hamiltonian (H) 
in the following ways: 

1. Holes cannot hop onto Ga. Therefore the average 
number of one third of a hole per unit cell, actually 
means that the fraction of Mn'^^ is increased with 
respect to the case ?/ = by a factor {1 — y)~^ . 

2. The core spins vanish at the Ga sites. 

3. A electrostatic potential appears on the Mn sites, 
proportional to the number of Ga'^"'" that are 
nearest-neighbors of the considered Mn site. 



III. RESULTS 

We have-studied our disordered model using the VMF 
techniqueJHj which is approximate for the spins, but 
treat charge carriers on each spin texture exactly (up to 
controlled numerical errors). The Hybrid Monte Garlo 
studyn3 has shown that the VMF overestimates the crit- 
ical temperature by a 30%, which is the same overes- 
timation factor found in classical statistical mechanics 
three-dimensional models with short range interactions. 
In the particular case of a disordered model, the VMF 
technique has the important merit of allowing the study 
of very large clusters (96 x 96 x 96 in this work). The self- 
averaging nature of the electronic density of states makes 
thus unnecessary to average our results over disorder re- 
alizations, which instead would be mandatory in a Monte 
Garlo study where the study of a single 16 x 16 x 16 clus- 
ter is at the, very limit of present day computers and 
algorithms.^ In order to make the study as unbiased as 
possible, we shall consider as VMF ansatzs all the spin 
ordcrings that were found in theJVIonte Garlo study of 
model in the absence of disorder ,Eil namely FM, G-AFM, 
A- AFM, C-AFM, skyrmion (SK), ffux, twisted and island 
phases. 

To gain some insight on the physics of our model, it is 
useful ±Q|Xpnsider (see Fig. |l|) the hole density of states 
(DOS)EIeI for the fully FM spin-polarized system with a 
20% Ga substitution, upon variation of the electrostatic 
parameter e. For e = 0, we have the usual DOS of the 
cubic lattice, slightly contracted and smoothed by the 
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FIG. 1: Hole density of states for a full spin-polarized system 
with a 20% Gallium (y = 0.2), and several values of the pa- 
rameter e. The curves are normalized to a maximum of 1 — j/ 
holes per unit cell, and the Fermi level is obtained integrating 
the DOS until reaching one third of a hole per unit cell. 



percolative disorder (holes cannot jump onto a Ga). We 
have a sharp contribution at midband, due to the local- 
ized states produced by the rare configurations of Mn 
sites fully surrounded by Ga sites (thereafter this very 
rare configurations were taken out by hand). If the elec- 
trostatic potential at a given site is large enough (in a Mn 
site with k neighboring Ga, it is fee), a hole state of en- 
ergy — fce will get localized in this Mn site, as the reader 
can check in Fig. |l|. This effect can he,investigated within 
the small defect- density expansion,E2l where e is allowed 
to be arbitrarily large. It turns out that a localized state 
will form if fee is larger than the edge of the e = D0S.l3 
Therefore, Fig. |l| tells us that for e = 2t, a Mn site should 
have at least three neighboring Ga to trap a hole on it, 
while for e = 3t it is enough to have two neighboring Ga 
(Mn sites with two neighboring Ga sites are deceivingly 
common, even for modest values of y, because of the 
large combinatorial factor). The electrostatic potential 
tends also to localize the electronic states even near to 
the edges of the main part of the DOS, depicted in Fig. |l| 
(for y = 0.2 the Fermi level is very close to the edge of the 
main part of the DOS). The calculations show that the 
double exchange interaction is strongly suppressed when 
holes occupy these localized states, leading to a decrease 
in the value of the Tq . We find that the drop in Tq upon 
Ga doping is best reproduced using e it, and, con- 
sidering 0U|LjaaX|Ii£x_results for manganites with a perfect 
Mn lattice,y'BBy we take Jaf ~ 0.06t— 0.07t. On the 
other hand, for e < 2t disorder is not the main agent of 
the eventual decrease of Tq as y increases (as happens in 
the case TR = Fe). 

In Fig. ^ we show the phase diagram of model (|^), for 
1/3 hole concentration, in the cases of y = (top) and 
y = 0.1 (bottom) Mn — Ga substitution. The very narrow 
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FIG. 2: Magnetic phase-diagram for our model at x — 1/3, 
as function of Jaf, for two values of the Ga concentration, y. 
Top: y = 0. Bottom: y = 0.1. 



phase-separation regions at y = 0.1 are left unlabelled for 
clarity. In agreement with the expected role of disorder 
near a discontinuous transition, Ga substitution induces 
a wedge of PM phase, down to the lowest temperatures, 
between phases with different types of long range mag- 
netic order. In the 0.06i — O.OTt Jaf range the SK and 
A-AFM phases are virtually degenerate. Between the A 
and C phases we find tiny regions of phase-separation, 
2x2x2 island phase and twisted orderings. Notice 
that the VMF assumes that the spin structure is spatially 
homogeneous, overestimating the tendency towards first 
order transitions in the disordered system. Thus, the re- 
sults shown in Fig. || strongly suggest (see below) the 
existence of a continuous transition at all temperatures, 
in the Ga doped system. 

In figure ^ (bottom) we show the phase diagram of our 
model, for four values of Jaf, namely Jaf = 0.05^, which 
is a reasonable value for La2/3Sri/3Mni_jyGaj,03 at least 
for y — 0^3 Jaf = 0.06i that marks the change from 
second to first order character of the PM-FM transition, 
and Jaf = 0.064t — 0.07i which we believe is an appropri- 
ate range for La2/3Cai/3Mni_j^Gay03. There are several 
points to be noted. First, see how at y = 0, the different 
values of Jaf are enough to explain the 30% differences 
in Tc between La2/3Sri/3Mn03 and La2/3Cai/3Mn03. 
As we see, for La2/3Sri/3Mni_yGay03, Jaf ~ O.Obt, we 
have a more or less monotonous decreases of the Tq, the 
ground-state being FM in all the interesting y range. On 
the other hand, the scenario is completely different for 
La2/3Cai/3Mni_2^Gaj;03 (Jaf = 0.064t and O.OTt). In 
this case, the FM ordering disappears for yc w 0.1 — 0.2, 
where the PM phase is stable until the lowest temper- 
atures. For larger values of y a very complicated situa- 
tion is reached with different, almost degenerated phases. 
This region might look somehow glassy experimentally. 
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FIG. 3; Bottom: phase-diagram for several values of Jap, 
including Jaf = O.OSt (La2/3Sri/3Mni_yGay03) and Jaf = 
0.064f — 0.07t (La2/3Cai/3Mni_yGay03), as predicted by our 
model. Top: phase-diagram for electrostatic potential e = 2t. 



The ordering temperature in the large y region is about 
one order of magnitude smaller than for y = 0. More 
remarkably, for y larger than i/c, the spins are not com- 
pletely polarized even at zero temperature (technically, 
this means that the mean-field at low tcmpeMture is aT, 
where a is a large but not infinite constantEa), and, also 
at zero temperature, the magnetization in the FM phase 
decrease upon approaching j/c although it does not get 
below 80% polarization. This effect will be amplified by 
the neglected quantum fluctuations of the spins, espe- 
cially at very low temperatures. Thus, given the wedge of 
PM phase that reaches near zero temperature at y^ and 
our underestimation of quantum fluctuations, it seems 
fairly plausible that a quantum critical point will be found 
in La2/3Cai/3Mni_j,Gay03, for y = 0.1 — 0.2. In or- 
der to emphasize how crucial the electrostatic attraction 
of holes to the Ga sites is, let us show in Fig. ^ (top), 
the equivalent of Fig. ^ (bottom), for e — 2t. In this 
case, the system is ferromagnetic in all the explored range 
< y < 0.4, in agreement with the scarcity of localized 
states shown in Fig. [I]. 



IV. REMARKS AND PREDICTIONS 

In this work we have accomplished a systematic 
study of doped manganites begun in previous pa- 
pers and formulated a microscopically motivated disor- 



dered model to study the effects of Ga substitution in 
La2/3Cai/3Mni_2^Gaj^03. We have studied this model 
using the VMF technique in 96 x 96 x 96 clusters. A result 
of this work is the prediction that a QCP is expected to 
appear for j/c ~ 0.1 — 0.2, which is induced by Anderson 
localization of the electronic states below the Fermi level 
that suppresses the double exchange mechanism. Thus, 
we expect localization effects in the transport properties 
at low temperatuxes, where phonons are static (ie po- 
larons could existo, but they would play the same role 
of our e diagonal static disorder term). For larger val- 
ues of y, a different phase is reached which is difficult 
to analyze because many different phases are almost de- 
generate. Nevertheless one can conclude safely that the 
ordering temperature for 0.3 < y < 0.4 will be in the 
20 to 40 K range. It is worth noting that the main fea- 
tures of the predicted phase-diagram are in quantitative 
agreement with preliminary experimental work.t^ 

A completely different behavior is expected for 
La2/3Sri/3Mni_j,Gaj^03. In such case the localization 
phenomenon that can be inferred from Fig. ^ for e = 3t, 
besides to the results of Fig. || corresponding to Jaf — 
0.05t, strongly suggest that the fact that holes are local- 
ized does not necessarily imply that the long range mag- 
netism disappears; in fact, to maintain it, it is enough 
that effective ferromagnetic interactions between close 
spins remain. Hence, for La2/3Sri/3Mni_j,Gay03 we ex- 
pect a FMmetal-FMinsulator transition at low tempera- 
tures (see also Ref. pO[ ). 

In the Ga case, it is worthwhile to remark that in our 
model, which considers spins coupled to charge carriers, 
the localization transition determines the magnetic crit- 
ical exponents too. Experimentally, it is very promis- 
ing that one be able to study a localization transition 
through a magnetic order parameter. If there is really a 
QCP for y « 0.1 — 0.2,Ea there would be presumably only 
one characteristic length in the system, which would set 
the scale both of the magnetic and transport phenomena. 

It is therefore interesting the theoretical and experi- 
mental determination of the magnetic critical exponents 
associated to the magnetic transition at very low T and 
its universality class. 

Finally, it is worth noting that the model used in this 
work includes the simplest combination of interactions 
compatible with the experimental data for the mangan- 
ites, in the absence of Ga doping. The confirmation of 
the predictions proposed here for Ga doped materials can 
illustrate the usefulness of simple models in studying the 
phase diagram of these compounds. 
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